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Strong longitudinal color flux fields will be created in the initial stage of high-energy 
nuclear collisions. We investigate analytically time evolution of such boost-invariant color 
fields from Abelian-like initial conditions, and next examine stability of the boost-invariant 
configurations against rapidity dependent fiuctuations[3]. We find that the magnetic back- 
ground field has an instability induced by the lowest Landau level whose amplitude grows 
exponentially. For the electric background field there is no apparent instability although 
pair creations due to the Schwinger mechanism should be involved. 

I. INTRODUCTION 

In the high-energy hmit of nucleus-nucleus collisions, the relevant degrees of freedom of the 
incident nucleus are the small-x partons described in the framework of the Color Glass Condensate 
(CGC) 2] with the typical momentum scale, or saturation scale Qs{x). On the other hand, success 
of hydrodynamic models in the data analyses at RHIC seems to imply the formation of quark-gluon 
plasma (QGP) in local equilibrium in a very early stage (t < 1 fm/c) of the event. Then it must 
be a decisive challenge to understand the physical mechanism for such a rapid equilibration from 
the CGC initial condition. 

Plasma instabilities especially the Weibel type, are intensively investigated so far, where one 
analyzes coupled equations of motion for the gauge fields and the hard particles with anisotropic 
momentum distribution. Despite that many intriguing phenomena have been reported, the kinetic 
description for particles with momentum p_L ~ Qa should be applicable only after particle formation 
T In the very early stage, t < 1/Qs, the system is so dense that the field description may 

be more appropriate. Here we shall investigate the instability problem in the classical Yang-Mills 
(YM) equations. 

The dense pre-equilibrium system appearing between the first impact and the equilibrated QGP 
was recently named Glasma^. The glasma, produced from the CGC initial condition, is still well- 
described by strong coherent YM field and can be treated as a weakly coupled system. Reflecting 
the CGC initial condition, this coherent field is boost-invariant in the high energy limit, and 
thus the time- evolution towards a thermalized QGP should involve the process exciting rapidity- 
dependent degrees of freedom. If rapidity-dependent unstable modes exist in the glasma, they will 
drive the system into an isotropic/thermalized state more efficiently. In fact, numerical simulations 
for the glasma have shown instability against the rapidity-dependent perturbations |5|. Here we 
present an analytic attempt for understanding the dynamics of the glasma instability [ij. 

II. EXPANDING COLOR FLUX TUBE 

In the CGC picture for the nuclear collisions, the random color sources are set on the light-cone, 
= (i ± z)/\/2, and cross each other at the origin (Fig. 1). The two sources are accompanied 
with pure gauge fields a\ 2(3^^) = x,y) behind, which uniquely fix the initial condition for the 

gauge field in the forward light-cone as (r = V2a;+a;~ and r] = \ ln(a;+/a;^))[6] 

-r^a{T - 0,x^) = r2|K(a;±), , (1) 

a*(r 0, x^) = a\{x±_) + a\{x±_) . (2) 




FIG. 1; Event setup (left) and schematic picture (right) of the color electric and magnetic flux tubes with 
transverse size 1/Qs, created between the two Lorentz-contracted nuclei just after the collision. 



Here we use the Fock-Schwinger gauge, At = 0. The corresponding initial field strengths are 

E'\r=o+ = ~ig[a{,a^] , ^1^=0+ = m^jK^ai] ■ (3) 

with all transverse components vanishing. Note that this setup is perfectly longitudinal-boost 
invariant and has no rj dependence while the transverse coherence length of the fields should be 
of order 1/Qs, reflecting the CGC structures of the incident nuclei. Thus we have the picture of 
multiple tubes of the boost-invariant color flux of size 1/Qs distributed in the transverse plane. 

One can compute the gluon production, in principle, by solving the source-free YM equations in 
the forward light-cone from this initial condition. Event averaging corresponds to averaging over 
the initial color source distribution in the CGC framework. Note, however, that we don't take this 
averaging, but presume that thermalization should occur in each event. In particular, we study 
here the time-evolution of a single isolated magnetic/electric color flux tube. 

Magnetic flux tube : A magnetic flux tube is the unique object in the glasma configuration. This 
can be realized with 

a(r, a:;j_) = 0, ai{T, xj_) . (4) 

For simplicity, we assume that the color direction of at and are the same and constant, and 
then one can ignore the commutators. The Gauss constraint is trivially satisfied in this case, and 
the YM equations reduce to the Bessel equation[l[. For definiteness, we take the Gaussian profile 
in the transverse plane for the initial condition: 

arHk^)o,~i^e~^ . (5) 

The solution is written in terms of the modified Bessel functions Iq^i{z): 

B%T>0,r) = Boe-^>\-Q"^"lo{2QlrT), (6) 
E^{T>0,r) = Bae-Q>\-Q>' h{2QlrT). (7) 

where = y/ (E^)"^ + (Ey)'^ and r = \x±\. This describes the outward expansion of the fiux 
tube as seen in Fig. 2 (left and middle). In the right panel of Fig. 2, we show the energy densities 
of the fiux tube integrated over the transverse plane, (B^)^ and {E'^)^ as a function of time r. 
Remarkably, this time dependence is quite similar to the numerical results in (j] . 



FIG. 2: Spatial profile of the longitudinal magnetic field (left) and the transverse electric field 
(middle), and time dependence of the averaged field strengths, (B^)^ and (iJ^)^ (right) for a single flux 
tube. The spatial profiles are plotted at five different times QsT = (solid), 0.5 (longest dash), 1.0, 1.5, 
2.0 (shortest dash) with setting Bq = 1. The initial value of (_B^)2(r) is also set to 1. Notice that — 
at T — 0. These plots are true for the electric flux tube after exchanging E and B. 



Electric flux tube : In contrast, an electric flux tube is realized with 

a{T,x±)^0, ai{T,x±):^0. (8) 

We assume again that the color orientations of a and i?^ are the same and constant. Then the YM 
equation reduces to the Bessel equation. With the Gaussian initial condition the following results 
are completely dual to the magnetic case by exchanging E and B. 



III. INSTABILITIES OF THE GLASMA 



The outward expansion of the boost-invariant flux tubes never achieves the isotropization of the 
system, because its longitudinal pressure is never positive (cx {B^Y — {E'^Y for the magnetic case). 
This is where the rapidity-dependent fluctuations play the leading part. 

We introduce the small fluctuations a;^^ as 

Ai= Ai{T,x±_) + a.i{T,ri,XA_) , A, = A,(t, + a,,(T, 77, a;_L) , (9) 

where Ai and Ar; are boost-invariant background flclds in the previous section. We perform stability 
analysis of this system against the perturbations ai^^[7.]. Note that coupling between Ai^ri and a^^^ 
is uniquely fixed in the non-Abelian gauge theory. For simplicity, we replace the background fields 
by r-independent and spatially constant electric/magnetic fields, and consider SU(2) group. 

Constant magnetic field : Setting the constant magnetic field B^ to the third color direction, we 
use the gauge, a = and = 6'^^{B /2){y6ii — x5i2)- By inspection, pure fluctuation is found 
stable, and now we study the pure fluctuations. After some algebra with help of the Gauss law 
constraint, one can derive for with the third color charge (±) and positive spin -I- [J, 0] 

ia.(r5.af )) + (^^ T mgB ± 2gB^ ~a^^^ + (^-dl + = , (10) 

where m and v, respectively, are the orbital angular momentum and the momentum conjugate 
to the rapidity (and a similar equation for negative spin — ). Note that the term ±2gB stems 
from the anomalous magnetic moment. Replacing the a:;_L-dependent part with the 2D harmonic 
oscillator wavefunction {Landau levels), the latter two terms in Eq. (10) are combined and yields 
(^ — gB)a"^^ for the lowest mode. At ^JgB r^ait = this "spring constant" becomes negative: 




FIG. 3: Re[7i„(2;)] for ^=8 (dashed) and 12 (solid), normalized at z=0.1. The lowest mode of Eq. (10) is 
given with the modified Bessel function U^l^/gB t). 

Nielsen- Olesen instability^,^. The typical scale ^/ gB ~ Qg. Interestingly, this means that the 
mode with v is stable until r^^ait — v/Qs, and then it grows up within Tgi-ow = Conversely, 
the maximum value i/max for the unstable modes at r satisfies QsT = Vmax + f ~ t'max- This may 
be related to the linear r-dependence of Vraax found in Q . 

Constant electric field : Setting the constant electric field to the third color direction, we 
take the gauge, a" = — (l/2)£'5°^ and ai = 0. In the case of a,, 7^ and Oi = 0, is found stable. 
Let us consider the opposite case: = and Ui ^ . Then the YM equations reduce to 

ia.(ra.6W) + jfci + ;! (^.± = 0, (11) 

where 6(t, k±) is related to the transverse part of di. One sees that the spring constant is non- 
negative and stable: no exponentially growing mode. From the physics point of view, we expect, 
and are to clarify, a certain relation of our formulation to the Schwinger mechanism for the particle 
pair creation. 

IV. SUMMARY 

We have studied the time-evolution of the boost-invariant color flux tube in a simplified Abelian- 
like setup, and found the outward expansion of the tube whose energy density evolves similarly 
to the numerical simulations ^4] . In the stability analysis we have found the unstable mode on 
the magnetic background but no apparent unstable mode in the electric case[l|, Q. To be more 
realistic, we need to treat the time-dependent background field, interactions among multi flux 
tubes, back-reactions of the fluctuations on the background fleld, and so on. 
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